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ABSTRACT 



Aims. A novel method of using hard X-rays as a diagnostic for chromospheric density and magnetic structures is developed to infer 
sub-arcsecond vertical variation of magnetic flux tube size and neutral gas density. 

Methods. Using Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) X-ray data and the newly developed X-ray 
visibilities forward fitting technique we find the FWHM and centroid positions of hard X-ray sources with sub-arcsecond resolution 
(~ 0.2") for a solar limb flare. We show that the height variations of the chromospheric density and the magnetic flux densities can be 
found with unprecedented vertical resolution of ~ 150 km by mapping 18-250 keV X-ray emission of energetic electrons propagating 
in the loop at chromospheric heights of 400-1500 km. 

Results. Our observations suggest that the density of the neutral gas is in good agreement with hydrostatic models with a scale height 
of around 140 ± 30 km. FWHM sizes of the X-ray sources decrease with energy suggesting the expansion (fanning out) of magnetic 
flux tube in the chromosphere with height. The magnetic scale height B(z) (dB/dzY 1 is found to be of the order of 300 km and strong 
horizontal magnetic field is associated with noticeable flux tube expansion at a height of ~ 900 km. 
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1. Introduction 

Chromospheric magnetic fields are notoriously difficult to mea- 
sure and their detailed structure is effectively inaccessible with 
modern observations. The difficulties of v arious line spectro- 
scopic techniques (e.g. ISolanki et al.ll2006l) occur because the 
magnetic field is relatively weak so the observed spectral lines 
are consequently broad and insensitive to the field. The compu- 
tation of chromospheric vector magnetic field fro m spectral lines 
is als o an ill-conditioned inverse problem (e.g. iMetcalf et alj 
I1995I) . In addition, current ground based vector magnetograms 
have spatial resolution compa rable with th e vertical size of the 
chromosphere itself 2 - 3" (Lagg 2007). Therefore, various 
indirect techniques are often employed to determine the mag- 
netic field in the chromosphere: optical observations of pho- 
tospheric magnetic fields combined with ex trapolation into the 
chromosphere (e.g. iMcClvmont et alj 1997b; radio observations 
of gyroresonance emission ( lLang et alJll993HAschwanden et ail 
ll995MVourlidas et al.|[l997HWhite & Kundull 19971) . 

The solar chromosphere being only about 2000 km thick 
(~ 3") strategically covers the layer where the solar atmosphere 
turns from the gas-dominated lower chromosphere/photosphere 
into the magne tic field dominated upper chromosphere/corona. 
iGabriell d 1976b has proposed th at magnetic field in the chr omo- 
sphere fans out (canopies) and iGiovanelli & Jonesl d 1982b have 
found that the canopy height should be t ypically 300-400 km. 
Howev er, polarisation measurements by lLandi deel'Innocenti! 
dl998l) suggest a very small hori z ontal component of the mag- 
netic field and ISchrijver & T itle (2003) argue that the "wine- 



glass" shaped magnetic field should return to the photosphere 
near their parent flux tube. In addition, different magnetic fi eld 
models predict different canopy heights (So lanki et al.lll999l) . 

The transport of both thermal and energetic charged particles 
in the solar atmosphere is governed by individual magnetic flux 
tubes. Therefore flare accelerated electrons one-dimensionally 
propagating along magnetic field lines can trace individual flux 
tubes from the electron acceleration site in the corona down to 
the deep layers of chromosphere where electrons emit hard X- 
ray emission. The asymmetry of hard X-ray footpoint sizes of 
a flaring loop (Melrose & White! 1 19791) could be a measure of 
the ratio of magnetic field strengths in the X-ray loop footpoints 
(ISchmahl et al.ll2007D . The simple dependency of X-ray emis- 
sion maximum location on the photon energy and density struc- 
ture in the low coro na/chromos phere due to Coulomb coll isions 
dBrown et al.l2002l) . has allowed lAschwanden etail d2002l) to in- 
fer the chromospheric density structure from Reuve n Ramaty 
High Energy Solar Spectroscopic Imager RHESSI dLin et al.l 
2002) X-ray observations. 

In this letter, we show that X-rays can be a diagnostic tool 
for the analysis of not only energetic electrons in solar flares but 
of the magnetic flux tubes and density structure in the chromo- 
sphere. We analyse the spatial and energy distribution of hard X- 
ray sources using RHESSI data to infer chromospheric density 
and magnetic field structure. Our results show that the density 
distribution of neutral hydrogen in the flaring loop has a scale 
height of 140 km and the magnetic flux tube of the flaring loop 
fans out by a factor of ~ 3 at the height of around 900 km. 
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Fig. 1. Cartoon of a magnetic loop and electron precipitation 
along the loop. The spatial size of the X-ray emission decreases 
with height, indicating the changing magnetic loop width. 

2. Electron precipitation in magnetic loops 

Accelerated electrons follow magnetic field lines towards the 
denser layers of the atmosphere, losing energy to binary colli- 
sions and scattering in an gle en route. In the first instance we 
follow iBrown et al.l (l2002l) and neglect scattering. Then the flux 
F(E,z) [electrons s _1 cm _2 keV _1 ] of electrons of energy E at 
depth z from the acceleration point z — 0, is 

F(E, z) = F in e( V£ 2 + 2KM3) 5 ~ \ ( 1 ) 

where N(z) is the column depth N(z) = J Q n{z)dz, K = 2n e 4 In A 
e is the electron charge, In A is the Coulomb logarithm dBrownl 
ll973HKontaretal.ll2002l) and we have assumed an injected flux 
of accelerated electrons F(E,z = 0) = F IN E~ S [electrons keV" 1 
cirT 2 s~'] at z — 0. The hard X-ray flux from depth z is 

/(e ' z)= iSIX F{E > z)oie > E}dE (2) 

where e is the photon en ergy and cr (e, E) is isotropic 
bremsstrahlung cross-section dHaud 11991) . Since the density 
n(z) increases with z while the electron flux F(E, z) decreases, 
I(e, z) thus exhibits a maximum for a fixed e at a depth z mflx (e) 
that increases monotonically with e (IBrown et al. I l2002h . Lower 
energy electrons lose their energy faster, with higher energy elec- 
trons propagating deeper into the atmosphere and the higher en- 
ergy X-ray emission should come from lower layers of the solar 
atmosphere. 

If the magnetic field strength B(z) increases with depth 
(Figure[TJ, the cross-sectional area A(z) of the flux tube will obey 
the principle of magnetic flu x conservation B(z)A(z) = const 
(e.g. iMelrose & WMtdfl979l) . The size of the X-ray source at 
energy e becomes a measure of A(z) at z ma x(e) fo r limb events 
(Figure [TJ, which are unaffected by albedo effects dKontar et al.l 
2006). We assume in the first instance that the electrons move 
parallel to the field so they experience no mirror force al- 
though B(z) varies. Electrons contributing most to radiation 
at Zmax(e) will still be close to their initial ener gies and thus 
mostly till moving along their initial directions (|Brownlll972l: 
iLeach & Petrosiadl 198 it MacKinnon & Craigll 19911) . 

3. RHESSI X-ray observations 

We selected a large GOES M6 class X-ray flare that appeared 
at the eastern limb on January 06, 2004 with a hard X-ray peak 
at ~06:23UT. The flare has an extended coronal source visible 
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Fig. 2. RHESSI CLEAN dHurford et alJl2002h X-ray images of 
the January 6, 2004 limb event. The contours show hard X- 
ray emission integrated for the impulsive phase of the flare 
(06:22:20-06:23:00 UT) from the footpoints in 30-35 keV 
(solid blue line) and 80-120 keV (dot-dashed green lines). The 
background image shows subsequent softer thermal emission 
(06:24:00-06:24:40 UT) in 10-18 keV. 

in soft X-rays and two hard X-ray footpoints (Figure 2). The 
"southern" footpoint is 5-10 times stronger than the northern 
footpoint, which seems partially occulted at energies above ~ 120 
keV. This flare seems ideal for our analysis since it is a limb 
event with one dominant source of hard X-ray emission seen 
in images up to ~300 keV. In addition, the spatially integrated 
count rate suggests emission above 800 keV. The spatially in- 
tegrated photon spectrum for the same time interval (06:22:20- 
06:23:00 UT) has been fitted using isothermal plus thick-target 
model (Figure 3). The X-ray emission above 18 keV (Figure 3) 
is dominated by the footpoint thick-target emission and can be 
used for our analysis. At these energies the "southern" X-ray 
source is the brightest and is used for detailed imaging analysis 
in different energy ranges. 

X-ray visibilities (2-dimen sional spatial Fourier compo- 
nents) (see ISchmahl et al.|[200~7b in ten different hard X-ray en- 
ergy ranges from 18 to 250 keV have been forward fitted us- 
ing single circular Gaussian source (Figure 3). Although the 
figure shows the fitted source matching the imaged footpoint 
location, the fit has been conducted on the X-ray visibilities 
from which the images are derived, not the images themselves. 
The X-ray i mages per se are generally a poor indic ator of the 
source size dHurford et al.ll2002t lEmslie et alj|2003l) . Visibility 
forward fit however provides reliable spatial measures and clear 
statistical uncertainties for all the fit parameters. These fits pro- 
vided us with the centroid positions of the Gaussian source 
(x(e), y{e)), and its Full Width Half-Maximum (FWHM) size 
s(e). Assuming vertically emerging magnetic field and calculat- 
ing the radial distance measured from the disk centre of the Sun 
R(e) = ^x(e) 2 + y(e) 2 we can readily find the height of the X- 
ray source h(e) = R(e)-Ro, where Rq is radial distance of the bot- 
tom of the loop. The typical uncertainties on radial distance mea- 
surement R(e) are around 0.2" ~ 150 km. Rq is poorly known 
(lAschwanden et all |2002|) but crucial for the detailed analysis, 
therefore we incorporate it as a fit parameter. Assuming the neu- 
tral hydrogen density profile to be 

n(h(e)) = no exp (~(R(e) - R )/h ) (3) 
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Fig. 3. RHESSI MEM NJIT (ISchmahl et all 120071) images in 
3 energy ranges (top to bottom: 20-22 keV, 35-43 keV, 55-75 
keV) integrated for the impulsive phase of the flare (06:22:20- 
06:23:00 UT). The FWHM circular source size found by for- 
ward fitting the X-ray visibilities, which were also used to make 
the background images, are shown by the overplotted contours 
(orange lines). 
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Fig. 4. Top: Spatially integrated RHESSI X-ray spectrum (his- 
togram) of the flare's impulsive phase (06:22:20-06:23:00 UT) 
fitted with an iso-thermal plus thick target model fit (solid line). 
The fit provides a temperature 1.65 keV (19MK), emission mea- 



sure 2 x 10 cm (dotted lines) and spectral indices 6 = 3.2 
below 191 keV and 3.9 above. Middle: "Southern" hard X-ray 
source height as a function of energy, h(e) = R(e) - Rq, where 
Ro = 975.3"; Bottom: FWHM circular source size as a function 
of energy. 



Figure 4) we find R = 975.3+0.2" and h = 140 + 30 km. Using 
these fitted parameters we can now plot the centroid height h 
against energy (Figure 4, middle panel), with the fitted model 
overplotted, finding that the height decreases by less than 500 km 
between 20 keV to ~ 200 keV. 

The FWHM sizes s(e) of the source, found from forward 
fitting the visibilities, decrease with energy from 6.2" (4.5 Mm) 
down to 2.3" (1.7 Mm) (Figure 4, bottom panel). Estimation of 
the magnetic field structure and density structure is shown in 
the top and bottom panels of Figure 5 respectively. We find that 
the magnetic field structure widens with increasing height and 
that the density decreases by over 2 orders of magnitude from a 
height of 1 Mm above the photosphere. 



where ho is the density scale height, and no is the number den- 
sity at height h - 0. Various chromospheric models A-F from 
IVernazza et al.l d!981l) provide slightly different values of hq 
though most are very close to no = 1.16 x 10 17 ctrT 3 . Therefore 
we use this no as a boundary condition to find the characteristic 
scale-height as well as Rq. Using equation ([3]) to fit the maximum 
of flux spectrum given by Equatio n (f2]i with In A = In A e H = 7.1 
jBrownlfl973t iKontar et alj|2002l) and using the spectral index 
found from the spatially integrated spectrum 6 — 3.2 (top panel 



4. Discussion and Conclusions 

Forward fitting X-ray visibilities allow simple and reliable mea- 
surements of not only locations of the emission maxima but 
also the characteristic sizes of hard X-ray sources. By using this 
visibility analysis on a good candidate limb flare we have im- 
pr oved earlier chromospheri c height and density measurements 
of lAschwanden et all (|2002|) . reducing the uncertainty of emis- 
sion maximum positions down to ~ 0.2". Assuming collisional 
transport in neutral hydrogen it can be concluded that the chro- 
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Fig. 5. Top: Chromospheric magnetic flux tube FWHM at var- 
ious heights with FWHM for density fit given by Equation [3] 
(solid line) Bottom: The chromospheric neutral number hydro- 
gen density as a function of height; best fit with hydrogen density 
profile given by Equa tion Q (solid line), density model C from 
IVernazza et all dl98ll) (dashed line). 



mospheric density is consistent with a gravitationally stratified 
atmosphere of density scale height 140 + 30 km. We show for 
the first time that not only is the higher energy X-ray emission 
produced continuously deeper in the chromosphere but the X-ray 
source sizes decrease from ~ 6.2" down to ~ 2.3". The precise 
measurement of the characteristic sizes allows us to conclude 
that magnetic field directing the energetic electrons converge, 
with the magnetic flux tube shrinking from FWHM ~ 3.5 Mm 
at h ~ 1 Mm down to FWHM ~ 2.5 Mm at h ~ 0.8 Mm 
above the photosphere. The magnetic scale height estimated us- 
ing flux conservation B(li) ~ FWHM(li)~ 2 as the ratio of areas 
BQi) (dB/dhT l = -FWHM(h)dh/dFWHM(h)/2 is found to be 
of the order of 300 km. This confirms the fanning ou t (canopies) 
of the chromospheric magnetic fields dGabrielll 19761) . 

Our deduced density structure agrees quite well with inde- 
pendent estimations from other methods, possibly a surprising 
conclusion in view of the simple treatment of electron trans- 
port, completely neglecting any pitch-angle changes. As already 
noted, collisional scattering will modify the range only by a 
factor of order unity but magnetic moment conservation might 
have a greater effect, unless fast electrons all start with veloc- 
ity vectors parallel to B. Large pitch angles would develop via 
collisional scattering only as electrons reached the end of their 
ranges, and the mirror force would be similarly unimportant 
most of the time. Such a concentration at small pitch angles 
seems at odds with findings of a nearly isotr opic electron dis- 
tribut ion from studies of photospheric albedo dKontar & Brownl 
120061) . The magnetic field convergence we find here however 
offers a simple resolution. The flux tube implied by our HXR 
source FWHM's implies magnetic field lines substantially in- 
clined to the vertical, possibly by as much as 60% on average at 
the outer edge and in addition are likely to be twisted. The HXR 
polar diagram of electrons populating the whole of an ensemble 
of field lines, with such a range of angles to the vertical, would 



be much closer to isotropic than expected for vertical field lines. 
The presence of magnetic canopies, then, appears to be a criti- 
cal factor for interpreting HXR images and even spectra, but one 
that needs a much more substantial and detailed modelling effort 
than we attempt here. 

While the X-ray diagnostics of the chromospheric magnetic 
field is an extremely attractive new technique, we emphasise that 
the morphology of X-ray sources should be further scrutinised. 
We note that the shape of the hard X-ray sources is more ellipti- 
cal (Figure 3) rather than circular suggesting that the vertical ex- 
tend of X-ray sources is govern by the spatially varying magnetic 
field (convergence and a twist of the flux tube), through the con- 
servation of the electron magnetic moment. This complicates the 
interpretation of the HXR source sizes beyond the simple ideas 
embodied in ([U and (O but brings a new diagnostic potential. 
This requires a more complicated electron dynamic and forward 
fitting elliptical sources to the X-ray visibilities, which will be 
the subject of future work. 

The method of measuring magnetic field using X-rays could 
be viewed as unique tool of loop width measurements in the 
chromosphere comparable to using TRACE data to measure th e 
widths of flux tubes in the corona dWatko & Klimchukll2000l) . 
While it is often argued that the heights of line formation should 
not be assigned due to the fact th at very distinct layers of the at - 
mosphere could be sampled (e.g. lSanchez Almeida et al]|1996l) . 
X-rays are uniquely related to the magnetic field lines connect- 
ing the photosphere and the electron injection site in the corona 
and weakly sensitive (via Coulomb logarithm dependency) to 
the temperature variations in the chromosphere. Therefore hard 
X-ray emission is likely to be a valuable diagnostic for mapping 
chromospheric magnetic field and density structures. 
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